EPJ manuscript No. 

(will be inserted by the editor) 



Diquark condensation in dense adjoint matter 



o 
o 

(N 

Or 
<D 

^sD 
(N 



> 
On 
(N 
O 
On 
O 



Simon Hands 1 , Istvan Montvay 2 , Luigi Scorzato , and Jonivar Skullerud 2 

1 Department of Physics, University of Wales Swansea, Singleton Park, Swansea SA2 8PP, UK 

2 Theory Division, DESY, Notkestrafie 85, D-22603 Hamburg, Germany 



43 ■ 



X 



the date of receipt and acceptance should be inserted later 

Abstract. We study SU(2) lattice gauge theory at non-zero chemical potential with one staggered quark 
flavor in the adjoint representation. In this model the fermion determinant, although real, can be both 
positive and negative. We have performed numerical simulations using both hybrid Monte Carlo and two- 
step multibosonic algorithms, the latter being capable of exploring sectors with either determinant sign. 
We find that the positive determinant sector behaves like a two-flavor theory, with the chiral condensate 
rotating into a two-flavor diquark condensate for /i > m^-/2, implying a superfluid ground state. Good 
agreement is found with analytical predictions made using chiral perturbation theory. In the 'full' model 
there is no sign of either onset of baryon density or diquark condensation for the range of chemical 
potentials we have considered. The impact of the sign problem has prevented us from exploring the true 
onset transition and the mode of diquark condensation, if any, for this model. 



1 Introduction 

In recent years, significant progress has been made in un- 
derstanding the phase diagram of QCD at non-zero baryon 
chemical potential (for reviews, see [0,0] )• On the basis of 
numerous model calculations, it is now believed that the 
ground state of QCD at high density is characterised by 
a diquark condensate which spontaneously breaks gauge 
and/or baryon number symmetries p,^|,|5],p|. However, al- 
though the results appear to be qualitatively independent 
of the specific model and approximation employed, lit- 
tle can be said quantitatively due to the lack of a first- 
principles, nonperturbative method that can access the 
relevant regions of the phase diagram. 

Lattice QCD, which would be such a method, fails be- 
cause the Euclidean-space fermion determinant becomes 
complex when a chemical potential is included, so stan- 
dard algorithms cannot be applied. However, it is possible 
to study QCD-like theories where the fermion determi- 
nant remains real even at non-zero fi. These theories can 
be used as testbeds to examine the validity of the models 
used to study real QCD, as well as directly to improve 
our understanding of phenomena such as diquark conden- 
sation and phase transitions in dense matter. Examples 
of such theories are two-color QCD, QCD with adjoint 
quarks or at non-zero isospin chemical potential 0, and 
the Nambu-Jona-Lasinio model ||,[|. A further point of 
interest is that chiral perturbation theory (yPT) may also 
be applied to many of these theories pO[[l]] ,p^|. 

Two-color QCD with a baryon chemical potential has 
been an object of study for lattice theorists for many years 
[]l3|,[l4|[l||]; there have recently been a number of simula- 
tions with quarks in the fundamental representation Jlq, 



H, p]§[l|,|20|, . Here we will be studying two-color 
QCD with staggered fermions in the adjoint representa- 
tion. The symmetries and conjectured phase diagram of 
this model have been presented in a previous paper [ pjj ; 
here we summarise the main features: 

— For an odd number of staggered flavors N, the fermion 
determinant may be negative. This means that this 
model still has a sign problem, which may make simu- 
lations at large /i difficult. 

— At m = /z = the U(iV) x U(N) flavor symmetry is 
enhanced to a XJ(2N) symmetry which relates quarks 
to antiquarks. This symmetry is broken by the chi- 
ral condensate to Sp(2iV), leaving N(2N — 1) massless 
Goldstone modes, which become degenerate pseudo- 
Goldstone states for m ^ 0. 

— When N > 2 these pseudo- Goldstone states include 
gauge invariant scalar diquarks, which are thus degen- 
erate with the pion at fi = 0. These models can be 
studied for // ^ by xPT jllj, the main result being 
that for /j > 771^/2 the chiral condensate rotates into a 
diquark condensate (the two being related by a U(2A) 
rotation) while the baryon density increases from zero. 
For N = 2 the diquark operator which condenses is 



QQs = - [X 



p tr 



(1) 



where p, q = 1,2 are explicit flavor indices and e the 
antisymmetric tensor. 

The model with A = 1 is not expected to contain 
any diquark pseudo-Goldstones and is not accessible 
to xPT. We expect an onset transition as some fi a sa 
mb/n q > to^/2, where is the mass of the lightest 
baryon and n q its baryon charge. 



2 



Simon Hands et al.: Diquark condensation in dense adjoint matter 



— For N — 1 the operator ([!]) is forbidden by the Ex- 
clusion Principle; there is, however, a possibility of a 
gauge non-singlet, and hence color superconducting, 
diquark condensate at large chemical potential: 

QQlc = \ [x tr (xYx{x) + x{xY l X tr (x)} , (2) 

where t l is a group generator in the adjoint represen- 
tation. Since qq\ c also transforms in the adjoint rep- 
resentation of the gauge group, the difference between 
the sum of the Casimirs of the constituents and that 
of the composite is positive, and hence the interaction 
due to one-gluon exchange is attractive in this channel. 

The indication in Q (and the update in is that 
the positive determinant sector of the N = 1 model be- 
haves like the model with N — 2. Ignoring the determi- 
nant sign has the effect of introducing extra 'conjugate 
quark' degrees of freedom q c , which carry positive baryon 
charge but transform in the conjugate representation of 
the gauge group (26) . In three-color QCD this approxima- 
tion leads to unphysical light qq c states which distort the 
physics of \i ^ beyond recognition. In two-color QCD 
the effect is more subtle; generically the physical spec- 
trum contains qq c states unless expressly forbidden by 
the Exclusion Principle. For models with staggered lat- 
tice fermions this is the case for N = 1 adjoint flavor p4j . 
In principle, therefore, by simulating this model we may 
get information about two 'physical' models for the price 
of one; N — 1 by taking account of the determinant sign, 
and N = 2 by ignoring it. The present paper aims to 
strengthen the evidence for this scenario. As discussed in 
p4j , for N = 1 staggered adjoint quarks, corresponding to 
four physical flavors, the continuum limit is problematic 
since the model is not asymptotically free. Our primary 
interest remains in studying a strongly interacting model 
with the potential to show superconducting behaviour re- 
gardless of these issues. 

The remainder of the article is organised as follows. In 
section || we study the performance of the two-step multi- 
bosonic algorithm. In section |j| we present our results. In 
sec. :h2 we study the chiral condensate, fermion density, 
and pion mass and susceptibilities, and demonstrate the 
different physical behaviour of t he positive determinant 
sector and the full theory. In sec. 3J3 we study pure gauge 
observables and the effect of Pa uli b lock ing in the posi- 



tive determinant sector. S ees. 0, 3^ and 3J3 all build on 
the results of p4| . In sec. |3.4| we present new results for 
diquark condensation, in both superfluid (i.e. gauge sin- 
glet) and superconducting (i.e. gauge non-singlet) chan- 
nels. Our conclusions are presented in section m 



2 Algorithm results 



The two-step multibosonic (TSMB) algorithm (27JH con- 
tains 6 tunable parameters: the polynomial orders n\ and 
n20; and the number of scalar heatbath, scalar overrelax- 

1 The orders , must merely be chosen large enough that 
they give rise only to negligible errors. 



ation, gauge metropolis and noisy correction steps Ih, lo, 
Im,I c in each update cycle. We have not systematically 
explored the effect of varying the relative values of Ih and 
Io, but we do not believe that this would have a major 
impact. In practice, we have also set I c = 1 which leaves 
us with 4 parameters to tune. 

The number of multiplications by M'M, where M is 
the fermion matrix, in a TSMB cycle is given by 

7 

"-tsmb ^ -^ni{I H + I + hi) + (n 2 + n 3 )I c . (3) 

The parameters ri\ and n 2 are essentially given by the 
condition number for a typical configuration at the simu- 
lation point. n\ is tuned to give a reasonable acceptance 
rate, while n 2 is tuned to give reasonable reweighting fac- 
tors. What is 'reasonable' in the latter case depends on 
whether we are in the phase where the determinant is 
positive, or whether there is a significant proportion of 
negative determinant configurations. In the former case, 
it is desirable to choose n 2 so large that the reweighting 
factors are very close to 1. In the latter case, the sign of 
the determinant will multiply the reweighting factor in 
the reweighting step, so a sharp peak around 1 will also 
give a sharp peak around —I. Instead, it is preferable to 
have a flatter distribution, thus permitting the algorithm 
to change the determinant sign. This is illustrated in fig. [j], 
where we show the reweighting factors obtained on a 4 3 x 8 
lattice with (3 = 2.0, m = 0.1 at three values for the chem- 
ical potential: /i = 0.3, where the determinant is always 
positive; fx — 0.37, where we find that 22% of the config- 
urations have a negative determinant; and fi = 0.4, where 
nearly half the configurations have a negative determi- 
nant. In all cases, we find that n 2 ~ 10ni. The optimal 
value of n\ increases somewhat with fi (and the condition 
number). 

In fig. |we show the distribution of the lowest eigenval- 
ues of the hermitean fermion matrix used in our update 
procedure at (i = 0.3 and 0.4. We see that the typical 
lowest eigenvalue, and thereby also the condition num- 
ber, changes by more than two orders of magnitude be- 
tween these two points. Also included is the lower limit 
e of the polynomial approximation we have used at these 
two points. We find that at [i = 0.3 the polynomial ap- 
proximation is nearly always accurate, while at (i = 0.4 
this is no longer the case for a substantial proportion of 
our configurations. 

How the autocorrelation depends on the values of Ih 
and Im is not obvious. It turns out that the acceptance 
rates are virtually unchanged when Im is changed. This 
means that it should be preferable to perform a long series 
of metropolis updates before the noisy correction, since 
the configuration will on average change more during the 
update sequence. This becomes inefficient at the point 
where the time taken for the metropolis updates begins 
to dominate over the time required for the noisy correc- 
tion. 

To determine the autocorrelation and its dependence 
on the algorithm parameters more precisely, we have per- 
formed long runs on a 4 3 x 8 lattice with = 2.0, m = 0.1 
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m=0.1 
|1=0.30 




0.5 1 -1.5 -1 -0.5 0.5 1 1.5 -1.5 -1 -0.5 0.5 1 1.5 

Reweighting factors Reweighting factors Reweighting factors 

Fig. 1. Reweighting factors, for fi — 0.3 (left), fi = 0.37 (middle) and [i = 0.4 (right). As the proportion of negative determinant 
configurations increases, the preferred distribution of reweighting factors changes from being sharply peaked around 1 to nearly 
flat. The polynomial orders (m,n 2 ) are (48,500), (80,800) and (100,1000) respectively. 

the presence of some very slow modes, and our runs are 
not long enough to accurately determine the autocorrela- 
tion in this case. The situation is substantially improved 
for Im = 20, although there are still substantial uncer- 
tainties and possibly residual slow modes in the data for 
/i = 0.37, where we have only been able to obtain order- 
of-magnitude estimates. We have combined all 8 estimates 
in one number or range for each parameter value; those 
are given in table |l|. 

We have also studied the autocorrelation of fermion 
quantities such as the chiral condensate and the pion prop- 
agator at fi = 0.3. Here we found an autocorrelation time 
of T/^;^ ~ ~ 120 with Im = 4, at least an order of mag- 
nitude smaller than for the plaquette. This means that we 
can be confident that the long plaquette autocorrelation 
times will not give rise to any additional uncertainty in 
fermionic observables. 













* hi,:- X*. ' r : .- : 4' '•V-/v;r'>-.C -s 











Fig. 2. The lowest eigenvalues of the hermitean fermion ma- 
trix, for fj, = 0.3 (left) and fj, — 0.4 (right). The lines indicate 
the lower bound e of the polynomial approximation employed. 



at /i = 0.3 and /i = 0.37. The parameters are given in 
table pi A measure of how expensive one update cycle 
(sweep) is, is given by the number Nih r of sweeps in one 
4-hour job on the Cray T3E, where this simulation was 
performed. We can see that the ratios of these numbers 
do not quite match those derived from (^); this may be 
attributed to the time the TSMB algorithm spends on op- 
erations that have not been included in the estimate (^). 

In all cases we have updated two separate lattices, and 
the autocorrelations have been measured separately for 
the timelike and spacelike plaquette for each lattice. The 
autocorrelation time has been measured by two methods: 
firstly by a straightforward measurement of the autocorre- 
lation function, giving the integrated autocorrelation time 
r a ; and secondly, by determining the jackknifc variance of 
the plaquette with successively larger bin sizes, until a 
plateau is reached. The autocorrelation time is then ob- 
tained as Tj f» var(n max )/2 var(l) where var(n) is the jack- 
knife variance with bin size n. In figure |^ we show the 
autocorrelation functions for each of our four runs. We 
can see that for Im — 4 there is a very long tail, and the 
autocorrelation itself is poorly determined. This reflects 



3 Physics results 

3.1 Simulation parameters 

As in (pl| , we have used three different masses using HMC 
at (3 = 2.0 on a 4 3 x 8 system, exploring values of /i up to 
and including 0.8 for m = 0.1, [i = 0.7 for m = 0.05, and 
/i = 0.5 for m = 0.01 p5[ . In addition we have performed 
further high statistics runs at m = 0.1 in the region of the 
onset phase transition at fi ~ 0.35. This latter parameter 
space was also explored using the TSMB algorithm, per- 
mitting an assessment of the impact of the determinant 
sign fluctuations. The simulation parameters in this re- 
gion for both algorithms are given in table |^. Note that in 
order to maintain a reasonable acceptance rate extremely 
short HMC trajectory lengths are required for /i > /i . We 
also performed a n a nalysis of diquark measurables to be 
described in sec. 3.4. 



3.2 Standard fermionic observables 

First we update our results for the 'standard' fermionic 
obervables: chiral condensate (ipip}? fermion density n and 
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Table 1. Parameters for our autocorrelation runs. In all cases we have used Ih = 2; Io = 8; I c = 1. n swccp is the length of the 
runs in update cycles; N&hr is the number of such update cycles achieved in 4 hours on the Cray T3E. titsmb is the number 
of matrix multiplications in one update cycle according to (H). r is the autocorrelation time in update cycles, while tm is the 
autocorrelation time in matrix multiplications. 
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Fig. 3. Autocorrelations for our four runs. The solid lines are the autocorrelations of the spatial plaquette; the dotted lines 
are the timelike plaquette autocorrelations. All these runs were performed independently on two processors, so there are two 
autocorrelation curves for each quantity. 



pion mass m n . In 24], wc found that the HMC algorithm, 
which samples only positive determinant configurations, 
yields results that agree well with the %PT predictions 
flllfl for theories with diquark Goldstone modes: 



nm v o 
%m{i>ijj)o 



1 x<l 

'o 

m^o x<l 
2[i x>l 



x<l 

X>1 



(4) 
(5) 
(6) 



where m is the bare quark mass, x = and the 

subscript denotes the values at /i — 0. Here we have 
extended our simulation to higher values of /i. The re- 
sults are shown in figs. ^, |^ and ^ for (ijjip), n and m v 



respectively. We see that good qualitative agreement with 
lowest-order xPT, in particular for (ipip), continues up to 
x ~ 2. It is also worth noting at this stage that the pre- 
diction for the onset x a = 1, clearly supported by our 
data, is stable to next-to-leading order in xPT p9[ . For 
x > 2, however, there is a dispersion between the data 
for different values of m. This could be explained by the 
appearance of new baryonic states in the spectrum not 
described by xPT (eg. a spin-1 qq state, or perhaps even 
a qg fermion), which should start to populate the ground 
state once /i ~ m},/n q . Since we expect such states to have 
a mass nib of order the constituent quark mass, governed 
by the magnitude of (ypip)®, then for different m the new 
thresholds should manifest themselves at different values 
of the rescaled variable x. For x > 3 there is also an in- 
dication of nonanalytic behaviour, which may conceivably 
be a sign of a further phase transition. It is also interesting 
to note that there are no signs of any saturation effects 
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Table 2. Parameters for the simulations in the neighbourhood of the onset. iV c fg is the total number of configurations used. 
A'gauge is an estimate of the number of independent gauge configurations, based on the plaquette autocorrelation times. With 
respect to fermionic observables, we assume all configurations to be independent, p- is the fraction of configurations with a 
negative determinant, while (r) is the average reweighting factor, with the sign of the determinant included. iVtraj is the number 
of HMC trajectories, and r the average trajectory length. 
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Fig. 4. Rescaled chiral condensate vs. chemical potential to- Fig. 5. Rescaled baryon density vs. chemical potential together 
gether with the x?T predictions ([|). with the x?T predictions (jB|). 



(i.e., when n approaches its maximum allowed value of 3), 
even at our largest /i = 0.8. It is clear, however, that any 
such conclusions remain provisional until data from larger 
volumes is available. 

In [^4| we found indications that the onset transition 
at x w 1 may be delayed when the negative configurations 
are included using the TSMB algorithm. Here we have 
collected additional data with both algorithms from the 
region just above this transition point for m = 0.1 (see 
table ||). The results are shown in figs. [7] and || for the chi- 
ral condensate and fermion density respectively. We can 
clearly see that the numbers for the positive determinant 
sector of TSMB agree well with the HMC data, giving the 
same early onset transition, but including the negative de- 
terminant sector brings the total average back so that it is 
consistent with the vacuum value — and inconsistent with 
the value in the positive determinant sector. We conclude 
from this that, as already indicated, the positive deter- 
minant sector and the full theory represent two different 



'physical' models: the former describes a two-flavor theory 
with an onset transition at fi = while the latter, 

the true N = 1 theory, remains in the vacuum phase for 
fi > rrin/2, and presumably has a transition to a normal 
or superconducting phase at higher fi. At /i = 0.4 we see 
that the uncertainties in the averages of the full theory 
become large. This is an effect of the sign problem, which 
becomes serious at this point. Simulating at even larger 
chemical potential, in the hope of finding the real onset 
transition in this model, would be very demanding and 
beyond our current resources. 

To gain further insight into the nature of the dense 
phase it is instructive to consider the quark energy density, 
given by 

e P hys = {x D oX) ~ |(3 -m(V^)o)- (7) 

The results are shown in table [|, along with the expected 
energy density nm^o/2 of a system with quark number 
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Fig. 6. m-n vs. fx, for the three different quark masses. Also Fig. 7. The chiral condensate from TSMB and HMC simula- 
tions. 
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Table 3. Energy density as a function of fi, compared with 
that of a system of non-interacting pions, for simulations with 
m = 0.1. 



density n consisting of non-interacting diquark baryons of 
mass m^Q. The two sets of numbers are comparable, al- 
though to the accuracy we have obtained it is impossible 
to assess even the sign of the binding energy in such a 
picture. It is also possible to express things in physical 
units; since m n oc -y/m J24j, and observing that our value 
of m/m^Q ~ 0.136 is about 5 times the physical ratio (as- 
suming m u = 4 MeV, = 140 MeV), we conclude that 
our vri — 0.1 simulation describes a world with m ~ 90 
MeV, 171^ ~ 670 MeV. This corresponds to a lattice spac- 
ing a ~ 0.22 fm. At (i = 0.37 this gives a quark number 
density n ~ 3.6 fm ~ 3 with a corresponding energy density 
e ~ 1270 MeVfm7 3 considerably in excess of the nuclear 
matter values n ~ 0.48 fm7 3 e ~ 150 MeVfmT 3 Perhaps 
a more reasonable comparison, however, is the dimension- 
less ratio e/(mn), which is ~ 3.8 at fi = 0.37 but has a 
value closer to 80 in nuclear matter. This emphasises the 
large value of the quark mass in our simulations; a more 
"realistic" simulation would require a separate calibration, 
e.g. via the vector meson mass. 
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Fig. 8. The fermion density from TSMB and HMC simula- 
tions. 



In table |4J we report the susceptibility defined by the 
integrated pion correlator 



E<^x(0)xex(z)> = GoMG&i-pi) , (8) 



where e(x) = (— \^ x o+xi+xi+xi an d G is the fermion prop- 
agator, which is real for adjoint quarks. In fact, the tab- 
ulated results come from 'non-singlet' diagrams with con- 
nected quark lines only; disconnected contributions were 
consistent with zero within large statistical uncertainties 
(we also attempted to measure the scalar susceptibility 
but did not find a significant signal with our statistics). 
The results are consistent, to within an irrelevant normal- 
isation factor of 3, with the axial Ward identity mx-w = 
(ipip). Once again the HMC data and the positive deter- 
minant sector of TSMB are in close agreement, showing a 
systematic decrease of Xn m the dense phase. Reweight- 
ing to the correct ensemble, however, again removes the 



Simon Hands et al.: Diquark condensation in dense adjoint matter 



7 





\Xir/TSMB 






/-w \ 


U.UU 


40.8 (( I ) 


















0.30 


45.76(5) 






45.87(22) 


0.35 








45.36(25) 


0.36 


45.79(24) 


45.28(21) 


35.32(155) 


45.04(17) 


0.37 


46.49(35) 


44.80(22) 


35.02(78) 


44.18(21) 


0.38 


46.23(81) 


42.49(30) 


36.34(56) 


42.91(45) 


0.39 








40.63(32) 


0.40 


46.24(321) 


38.06(34) 


35.56(38) 


37.42(40) 


0.50 








21.15(82) 



Table 4. Pion susceptibility x-n from TSMB and HMC simu- 
lations. 
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Table 5. The pion mass m n from TSMB and HMC simula- 
tions. 



evidence for the onset transition - this is significant since 
it demonstrates that reweighting also gives the expected 
behaviour for two-point correlation functions. 

The pion mass m w (/Lt) is extracted from the temporal 
decay of the correlator over timeslices 1-7. The results are 
shown in table ||. HMC results show that m n is constant 
and equal to m^o for fj, < [i , and increases for fj, > \x a in 
approximate agreement with the x^T prediction ([|); the 
fact that the results lie slightly above the xPT value is 
probably attributable to the small volume. It is interesting 
to contrast this behaviour with that observed in two-color 
QCD with fundamental staggered quarks p2| , in which 
rriir is observed to decrease once [i > fi . This behaviour is 
predicted by %PT Q for a pseudo-Goldstone meson with 
a flavor content symmetric under the residual global sym- 
metry in the superfluid state, which for N — 2 staggered 
adjoint flavors is Sp(2); m v should decrease with [i in a 
theory with Dyson index j3 = 4 such as two-color QCD 
with staggered fundamental quarks, and increase if /3 = 1 
as is the case for the current model. It becomes very dif- 
ficult to determine from the TSMB simulation once 
/i > 0.3, so we are not able to draw any significant conclu- 
sions from these data. However, the general tendency is 
compatible with what we observe for (tpij)) and n; namely, 
that m-n- increases (with values compatible with those from 
HMC) when only positive determinant configurations are 
included, while it remains compatible with the /j, = value 
when all configurations are included. 



3.3 Gauge field quantities 

Our results for the plaquette and Wilson loops are given in 
table ^. Beyond the observation that both algorithms give 
consistent results for the plaquette we see no discernible 
effect over this range of chemical potentials — all quanti- 
ties remain consistent with their values at /i = 0. At larger 
fi there is evidence that the plaquette starts to fall towards 
its quenched value due to the effects of Pauli blocking 
p5| . The deviation of the /i = 0.37 values from the rest is 
probably a sign of insufficient statistics and/or a too short 
equilibration time. 

The Polyakov loop L on small lattices in the dense 
phase of two-color QCD with fundamental quarks was 
found to be small but non-zero in p^| . Our results for 
the average Polyakov loop are also given in table |]. We 
can see that it remains zero everywhere, and conclude that 
there is no sign of any deconfinement transition either in 
the positive determinant sector or in the full theory. 



3.4 Diquark condensation 

A natural explanation of the agreement between the HMC 
and positive determinant results and xPT predictions is 
that the positive determinant sector mimics a theory with 
two flavors: one quark and one 'conjugate quark' (g6|. We 
may cast further light on this by studying the behaviour of 
the diquark condensate ([!]) expected in the two-flavor the- 
ory. This condensate may be determined by introducing 
a diquark source term in the action, which now describes 
two flavors, 



S\j] = S + j^2qq 3 (x) . 



(9) 



extracting the condensate {qqa{j)) =V~ 1 dlnZ[j]/dj, and 
extrapolating the results to j = Jl6|. Here, we have 
only performed 'partially quenched' measurements — i.e., 
we have included a nonzero diquark source j in the mea- 
surement of observables from configurations generated at 
5=0. 

In figure |9| we show the chiral and two-flavor diquark 
condensates from our HMC simulations as a function of 
j for a number of values of the chemical potential. Note 
that on a finite system at j = the diquark condensate 
is identically zero. We have therefore excluded this point 
from our plot. 

The extrapolation j 1 — ► needs some discussion. To 
leading order in the relation between chiral and di- 
quark condensates, m, and j is pj 



x 2 (ipip}(qq 3 ) = 



together with the constraint 



m 



(10) 



(ii) 



If we assume j <C m and ignore quadratic terms, then at 
at the critical point x Q = 1 the resulting cubic equation 
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0.1377(39) 
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0.1478(30) 


0.1477(24) 


0.1474(25) 


0.40 
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L 

0.00 


-0.003(3) 






0.30 


-0.005(4) 






0.36 


0.015(9) 


0.014(9) 


-0.002(10) 


0.37 


-0.002(4) 


0.000(4) 


0.011(7) 


0.38 


0.010(9) 


0.008(7) 


0.004(10) 


0.40 


-0.005(18) 


-0.010(6) 


-0.11(6) 



Table 6. Results for the gauge observables: the plaquette □, 
spatial plaquette □ timelike Wilson loops Wij, and Polyakov 
loop L. 
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Fig. 9. The chiral and diquark condensates as a function of 
the chemical potential fi and the diquark source j, from HMC. 



yields (qq a )/(H>)o =W&)* - §£• 
tried to fit the data using the form 



We have therefore 



(qqati)) = A + Bp +C] . 



(12) 



The results of the extrapolation are shown in fig. |l0[ Pos- 
itive values for the coefficient A were found for fi > 0.36, 
and the quality of the fit improved as fi increased. At 
fi = 0.36 the values for the coefficients were A — 0.062(15), 
B = 9.37(7), C = -16.1(2), the latter to be compared 
with xPT predictions at x = 1 of 8.28 and -20.4 respec- 
tively. Although the quality of the fits improves as fi in- 
creases, the curvature due to the j's term is observed to 
manifest itself principally at values of j < 0.01, i.e., be- 
low where we have data. We therefore suspect that fits to 
( |l2|) underestimate (qqsU = 0)) for [i > [l . This is in line 
with the solutions of ( |lQ|Jll| ) at a; > 1: extrapolating the 
solutions of (|l0|) at j > to j = using ([l2]) always yields 
a lower value than the actual solution at j = 0, 



(qqs) = WOoyi 



1 



(13) 



We can also see from fig. ||that the behaviours of (ipip) 
and {qqz) are strongly anti-correlated. This is a result of 
the two quantities being connected by a U(4)/ transfor- 
mation 124(1, 



X = 



Xo 

xt 



VX X = (XeXt 



XV 



(14) 



Simon Hands et al.: Diquark condensation in dense adjoint matter 



9 



1 1 1 1 1 1 1 


I 


= ° ° T 

/"V 


• 

/' A 


\ 

o<w 3 > Y 

A<qq 3 > MFT extrapolation ^ 




I i 1 i I A I 


i 



0.1 



0.2 



0.3 



0.4 



0.5 



Fig. 10. (V't/'), (993) and G versus fi, together with predictions 
from leading order x?T. "MFT extrapolation" denotes points 
from fits to (Til) 



with 



V 



piptr 
iP P tr 



where P 



(15) 



The superfluid transition in this model is expected to be 
realised by the chiral condensate 'rotating' into the di- 
quark condensate, in such a way that the constraint ( [TT| ) 
is maintained. This has been shown to be the case to lead- 
ing order in xPT |pj| . At next to leading order, however, 
some dependence on ji and j is expected j2{j] . 

We have therefore tried a second method in which the 
quantity G 2 — (ipip) 2 + {qqz) 2 measured for j > is 
extrapolated to j = using a second order polynomial. 
Fig. [n] shows that G^j, n) exhibits relatively little varia- 
tion with either j or [i, implying that the NLO xPT cor- 
rections are small. Note that since (qqsij = 0)) =0, the 
extrapolated G^O) lies considerably above the numerical 
data for fi > fi a . We now use G(fJ), which appears to vary 
very little across the transition at /i = /i , together with 
($V>) to extract (qqa); the results are plotted in fig. [l0| In 
this case difficulties associated with taking the difference 
of two fluctuating quantities means that the method is 
probably not accurate for {qqs) in the immediate vicinity 
of the critical point; for /i > 0.37, however, the agreement 
with the yPT predictions (p|,|l3|) is striking. 

In figure [l2| we compare the HMC results for {qqs) at 
/U = 0.38 with TSMB numbers. Since the superfluid tran- 
sition is not present for N = 1, the diquark condensate 
should be suppressed, and indeed the data from the full 
theory lie lower and exhibit a stronger negative curvature 
than the positive determinant data. It is, however, not pos- 
sible to conclude from these data whether the one set will 
extrapolate to zero and the other to a non-zero value, es- 



G 8 
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Fig. 11. The quantity G 2 (j) for values of p in either phase. 
Lines show the quadratic extrapolation to j — 0. 
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Fig. 12. The diquark condensate as a function of the diquark 
source j at fx = 0.38, from TSMB and HMC. 



pecially taking into account the problems discussed above. 

We have also computed the following four-point func- 
tions, which we refer to as superfluid and superconducting 
"susceptibilities" : 

X3 = {x{x)ex tr {x)x tr {x)e X {x)), (16) 

Xsc - l(x(x)fx tr (x)x tr (x)f X (x)) . (17) 

Both may be computed without introducing a diquark 
source term in the action |||l6), and are related to the 
superfluid and superconducting diquark condensates of 
eqs. (0) and respectively. Both are scalar objects cor- 
responding to the local component of a diquark suscep- 
tibility Xqq — J2x(l < l(^)'l'l( x )) w hich we might expect to 
increase significantly if condensation occurs. Note, how- 
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Fig. 13. The superconducting (top) and superfluid (bottom) diquark susceptibilites as a function of the chemical potential fi, 
from HMC and TSMB. 



ever, that Xsc in (^7|) is the only such contribution which 
is gauge invariant and hence a possible signal of color su- 
perconductivity in a simulation (to be compared with <P^<P 
in a Higgs model). To be able to compare signals for con- 
densation between channels we have therefore chosen the 
local operators ([l6 17) over the full susceptibilities. 

Figure [l3] shows the superfluid and superconducting 
diquark susceptibilities as a function of chemical poten- 
tial from our HMC and TSMB simulations. In the low 
density regime xa an( A Xsc show little variation with /x 
and are roughly equal. Once we enter the superfluid phase 
for /x > fi a , however, signals for both observables become 
increasingly dominated by very sharp peaks of amplitude 
0(10 — 20), the distinction being that for xa the peaks are 
all positive whereas for Xsc they occur with either sign. 
Such peaks are characteristic of localised small-eigenvalue 
modes of p, of which there are many in the dense phase 
p4| . For fi > /i , in the positive determinant sector Xsc 
decreases. We interpret this as an effect of Pauli blocking 
or phase space suppression: as the ground state is filled up 
by fermions, there is less phase space left to accommodate 
the fermion loops that contribute to Xsc- 

When the sign is taken into account, this effect, like 
other effects of the superfluid transition, disappears. Per- 



haps surprisingly, this is in fact the cleanest signal we have 
of the effect of the sign of the determinant on observables. 
By way of contrast, xa shows a small increase as the onset 
transition is reached, as one would expect as a result of a 
Goldstone mode once a superfluid condensate forms. How- 
ever, there is a large deviation between X3 as measured in 
HMC and TSMB(+); this is the only observable we have 
examined in which the two algorithms significantly fail to 
agree. One possible explanation is that Goldstone physics 
is extremely sensitive to small eigenmodes of lp. In HMC 
such modes are responsible for large force terms and re- 
duced acceptance rates; moreover the sampling may not 
be correct due to the breakdown of ergodicity . Hybrid 
algorithms have produced disparate results for susceptibil- 
ities in the dense phase |l(||l7]]. TSMB, in contrast, tends 
to over-sample configurations with small modes, correct- 
ing this by reweighting factors with modulus less than 1. 
In fact, in this case some of the eigenvalues and reweight- 
ing factors are very small, and therefore problems with 
numerical accuracy cannot be excluded. It is perhaps not 
surprising that the algorithms disagree in this exacting 
regime. We were not able to obtain any result for the 
overall superfluid susceptibility in TSMB (i.e., with the 



Simon Hands et al.: Diquark condensation in dense adjoint matter 



11 



negative determinant sector included) due to huge fluctu- 
ations in the signal. 



4 Conclusions 

We have studied two-color QCD with one flavor of stag- 
gered quark in the adjoint representation. We have em- 
ployed two different simulation algorithms, and have con- 
tinued to gain insight into the optimal tuning of the TSMB 
algorithm in the high density regime. This is t he p referred 
algorithm not only because as highlighted in [E4| it is ca- 
pable of maintaining ergodicity via its ability to change 
the determinant sign once )jl > /x a , but also because it 
more effectively samples small eigenmodes, which are im- 
portant in the presence of a physical Goldstone excitation. 
We find that the positive determinant sector behaves like a 
two-flavor model, and exhibits good agreement with chi- 
ral perturbation theory predictions for such a model in 
the regime 1 < x < 2. At higher chemical potentials there 
are preliminary indications of a breakdown of xPT, and 
possible signs of a further phase transition. However, data 
from larger volumes and smaller bare quark masses would 
be needed to make these observations definitive. 

Above the onset transition in the positive determinant 
sector, we have successfully obtained a signal for a non- 
zero two-flavor diquark condensate (qqs), indicating a su- 
perfluid ground state for [i > The chiral conden- 

sate rotates into this diquark condensate, in good quan- 
titative agreement with the behaviour predicted by xPT. 
This feature also enabled us to achieve reasonable con- 
trol over the necessary j ' — > extrapolation. We also find 
that the superfluid susceptibilit y ( f Lq ) increases, while the 
superconducting susceptibility (|17[) decreases for /i > fi Q . 
For the former there are indications of a lack of consis- 
tency between HMC and TSMB algorithms, indicative of 
the two methods' differing treatment of small eigenmodes. 
The latter may be interpreted as an effect of phase space 
suppression. Unfortunately we have seen no evidence for a 
superconducting condensate (qql c ) ^ 0, one of our original 
motivations for studying the model. 

When the negative determinant configurations are in- 
cluded in the measurement, the onset transition and di- 
quark condensation disappear. This is what we would ex- 
pect for the one-flavor model and is consistent with the 
Exclusion Principle. There is much stronger evidence for 
this scenario than that presented in our previous paper 
p4j , providing a conclusive demonstration, should any still 
be needed, that the determinant sign plays a decisive role 
in determining the ground state of systems with fi =/= 0. 
Unfortunately, the severity of the sign problem means we 
have not been able to locate the real onset transition for 
this model. 
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